Abstract. Three types of observations: the daily values 
of the solar radio flux at 7 frequencies, the daily interna- 
tional sunspot number and the daily Stanford mean solar 
magnetic field were processed in order to find all the pe- 
riodicities hidden in the data. Using a new approach to 
the radio data, two time series were obtained for each 
frequency examined, one more sensitive to spot magnetic 
fields, the other to large magnetic structures not connected 
with sunspots. Power spectrum analysis of the data was 
carried out separately for the minimum (540 days from 1 
. March 1996 to 22 August 1997) and for the rising phase 
^ ; (708 days from 23 August 1997 to 31 July 1999) of the so- 
fS| ■ lar cycle 23. The Scargle periodograms obtained, normal- 
^ ized for the effect of autocorrelation, show the majority of 
known periods and reveal a clear difference between the 
^0 periodicities found in the minimum and the rising phase. 
Qv^ . We determined the rotation rate of the 'active longitudes' 
^ ; in the rising phase as equal to 444.4 ± 4 nHz (26':'0±0'?3). 
The results indicate that appropriate and careful analysis 
of daily radio data at several frequencies allows the inves- 
tigation of solar periodicities generated in different layers 
' of the solar atmosphere by various phenomena related to 
, the periodic emergence of diverse magnetic structures. 

0\ ; 
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1. Introduction 

The study of periodicities in different solar data is impor- 
tant for the understanding of solar magnetic activity. For 
over a decade many authors have reported various periods 
other than those at 11 yr and 27 days, of which the one 
near 154 days is the best known (Rieger et al. 1984; Dennis 

Bai & 



1985; Bogart & Bai 1985 



Cliver 1990; Droge et al. 



& Ballester 1990, 1992 



Bouver 



Bai & Sturrock |1991 
Oliver et al 



1992 



Lean & Brueckner 1989 
1990| Pap et al. |199C| ; CarboneU 



199S; Ballester 



& Cli ver |1991| ; 

et al. 1999| ). This investigation was designed to show that 
accurate daily radio observations can be very useful for 
studying the solar periodicities generated in different lay- 
ers of the solar atmosphere, and in different phases of the 
solar cycle. We used a new approach to the daily measured 
radio fluxes which allow us to separate to some extent the 
radio emission generated in the strong magnetic fields of 
active regions from that emitted by large but weaker mag- 
netic structures. All the solar data from the minimum and 
the rising phase of solar cycle 23 were processed sep aratel y 
using the Scargle periodogram technique (Scargle 1982| ). 
The statistical significance levels of all the periods found 
were estimated through the FAP (false alarm probabil- 
ity - the probability that the given periodogram value z 
is generated by noise), as well as by a Monte Carlo ap- 
proach to the periodograms obtained (Horne & Baliunas 
1986|; Bai |l992b| Ozgiig & Atag |1994|; Oliver & Ballester 



1995) 



Papers published by Das & Chatterjee ( |1996D ; Das & 
Nag (1998, 1999) presented some periods in radio data 
without any discussion about their significance level. Our 
analysis of radio observations however, gives the majority 
of the known periods, reveals a clear difference between 
periodicities observed in two phases of the solar cycle and 
shows that daily measured radio fluxes at various frequen- 
cies are very useful for the systematic study of solar pe- 
riodicities observed in the different layers of the solar at- 
mosphere. 
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2. Observations 



The radio fluxes at the four lowest frequencies 405, 
810, 1215, 1620 MHz are from our new accurate 
observations, which we began on 1 March 1996 

Every day we measure the to- 



( http: / /www.oa.uj .edu.pl 



199^ ; Kile et al 



tal solar radio flux at ten frequencies within the decimeter 
range of wavelengths using the 8-m radiotelescope built 
in 1995 at the Cracow Astronomical Observatory (Zi§ba 
1996| ). The international sunspot numbers {ISN), 



the Stanford mean solar magnetic field (MMF) and radio 
fiuxes at the frequencies 2800, 4995 and 880 MHz wer e 
taken from Solar Geophysical Data bulletins ( 1996-1999 ). 
We added to our radio observations those at higher fre- 
quencies for two reasons. Firstly, we wanted to examine 
periodicities of the radio emission in a wide enough ra- 
dio band to analyse the radio flux coming from a large 
part of the solar atmosphere and secondly to compare our 
observations with those of others. 

It is trivial to form the time series from the daily val- 
ues of the sunspot numbers and the mean magnetic field, 
but if we want to use the radio data for the study of 
magnetic activity, we must first eliminate from the ob- 
served daily flux (free from bursts) the thermal emission 
which comprises the majority of the daily measured value. 
This can be done through a procedure in which the flux 
observed on day ti, Fi is divided between the thermal, 
almost constant component (often called the basic com- 
ponent - B) and the slowly varying component {SVC)^ 
whose value changes every day and is generated by mech- 
anisms dependent on the magnetic field. It is usually as- 
sumed that this component is proportional to a certain 
daily index of activity, for example {ISN)i and then the 
daily fiux Fi can be described by the following linear for- 
mula: Fi = B + {SVC)^ = B + h{ISN), , where h is the 
production of the radio flux from the spot with ISN = 1 
(Kruger p^79| ). 

The assumption that B is constant over the large time 
intervals is rather strong and cannot be accepted, espe- 
cially at a time when the level of radio flux rises system- 
atically. Here we propose a new approach in which the 
basic component B is not constant over the given time 
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Table 1. Parameters of the two different models of the solar radio emission used to eliminate the basic component B 
from the daily observed flux Fi. 



Linear model Fi = B + h{ISN)i 


frequency in MHz 


405 


810 


1215 


1620 


2800 


4995 


8800 


B- basic component [suj 


22.9 ±.2 


39.8 ±.3 


48.2 ±.4 


53.5 ±.4 


67.3 ±.5 


117.0 ±.5 


199.4 ±.1 


h- radio flux production 


0.138 ±.003 


0.337 ±.005 


0.508 ±.007 


0.611 ±.007 


0.785 ±.008 


0.633 ±.008 


0.752 ±.017 


% of the data variance 


71.0 


78.6 


81.9 


85.7 


87.5 


83.5 


61.6 


explained by the model 
















Boltzman formula model 


Fi = Bi + hs 


[ISN)i = A2 + 


(Al - A2)/(l + exp((t, - to)/ At)) + hsilSN), 






frequency in MHz 


405 


810 


1215 


1620 


2800 


4995 


8800 


parameters of Ai [suJ 


22.1 ±.5 


38.9 ±.8 


48.0 ±1 


53.1 ±1.2 


67.4 ± 2 


117.8 ± 1.4 


200.5 ± 7.6 


the Boltzman A2 [su] 


31.4 ±.5 


62.5 ±.7 


75.0 ±.9 


78.9 ± 1.1 


100.7 ± 1.8 


137.1 ± 1.5 


292.5 ± 7.1 


sigmoidal to [days] 


807 ±15 


855 ± 7 


879 ± 8 


880 ± 10 


926 ± 12 


962 ± 11 


1064 ± 20 


formula At [days] 


132 ± 13 


111 ± 7 


80 ± 7 


91 ± 9 


106 ± 12 


45 ± 12 


143 ± 16 


Hb- radio fiux production 


0.067 ±.004 


0.154 ±.006 


0.291 ±.008 


0.408 ±.010 


0.542 ±.012 


0.487 ±.012 


0.261 ±.025 


% of the data variance 


81.0 


91.8 


91.4 


91.5 


92.5 


87.5 


79.3 



explained by the model 



Cracow radio observations 
at 810 MHz 



Rising pliase 
23 Aug 1997 -31 Jul 1999 



Minimum 
1 Mar 1996 - 22 Aug 1997 




B = Boltzman formula 



Basic component B = constant 



150 300 450 600 750 900 1050 1200 1350 
Time (days alter 1 Jan 1996) 

Fig. 1. Daily values of the radio flux at 810 MHz observed 
from Cracow. The horizontal solid line shows the constant 
value of the basic component B = 39.8 su resulting from 
the linear formula, while the dashed curve shows values 
of the basic component calculated from the best fitted 
parameters, Ai = 38.9 su, A2 = 62.5 su, to — 855 days, 
At = 111 days according to the Boltzman formula Bi ~ 
^2 + (^1 - ^2)/(l + exp((i^ - to)/ (At)). The division into 
the minimum and rising phase is also indicated. 

interval but changes every day according to the Boltzman 
sigmoidal formula: 



B, 



A2 + (Al - yl2)/(l + exp((t, - to)/At)) 



where Ai, y4.2, to, At are the parameters and t^ is the day 
number. 

Then the observed daily flux Fi ^ Bi + h^{ISN)i, 
and /iB has a similar interpretation to h. To determine 
the above parameters we used the observed radio data 
and daily sunspot numbers over the whole time interval 
investigated, 1 March 1996 - 31 July 1999 (1248 days). The 
best fit values of these parameters are shown in Table 
1. The difference between the two models is clearly seen. 



especially for four frequencies 405, 810, 1215, 8800 MHz. 
To demonstrate this we present in Fig. 1, as an example, 
the daily values of the radio flux at 810 MHz observed 
from Cracow as well as the calculated values of the basic 
component B and Bi. The data in this flgure also explain 
our division into the minimum and the rising phase. 

Thus, in our approach to the radio data we create two 
time series from the observation at each frequency. The 
first, the SVC (slowly varying component) time series, 
consists of diurnal values calculated as the difference be- 
tween the daily observed ffux and the daily value of the 
basic component computed from our model, {SVC)i = 
Fi — Bi. The second, the RRE (radio residual emission) 
time series describes the every day difference between the 
radio observations and our model of the daily radio ffux, 
(RRE), = Fi - Bi - hB{ISN)i. Taking the time series 
SVC and RRE, we can analyse cyclic variations of those 
magnetic structures which modified the observed radio 
emission. However, the SVC series are more sensitive to 
spot magnetic fields, while the RRE series are sensitive to 
large magnetic structures not connected with sunspots. 



3. Periodogram analysis 

The search for periodicities in all the time series was 
performed by calculating the Scargle normalized peri- 
odograms PAr(w) (Scargle 1982| ).This technique (see Horne 
& Baliunas 1986| ) has several advantages over the conven- 
tional fast Fourier transformation method and provides, 
through FAP, a simple estimate of the significance of the 
height of a peak in the power spectrum. However, the 
FAP value is easy to calculate only for that time series 
for which the successive data are independent. In our case, 
the all analysed time series were prepared from the daily 
values of the different solar indices which are not inde- 
pendent but correlated with a cha racter istic correlation 
time of a week (Oliver & Ballester 1995 ). Therefore, for 
all our FAP calculations used to estimate the statistical 
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Scargle power log P, (Period P in days) 

Fig. 2. a). The cumulative distribution function of the Scargle power for the original, minimum ISN time series. The 
vertical axis is the number of frequencies whose power exceeds z. The straight line is the best fit to the points for 
values of power lower than 5. b). The normalized periodogram of the original, minimum ISN time scries with FAP 
significance levels indicated. 



significance of a peak z, in the Scargle power spectrum, 
we applied the formula FAP = 1 — [1 — exp(— z/fc)]^ = 
1 — [1 — exp(— Zm)]''^, where N is the number of indepen- 
dent frequencies, k is the normalization factor due to data 
correla tion and — z /k the normalized power (Bai & 
Oliver |l990| ; Bai |l992b|) . To determine the normalization 
factor k for the given time ser ies we followed the proce- 
dure described by Bai & Oliver ( 1990 ). The key step of this 
procedure is the choice of a spectral window as well as the 
number of independent frequencies. Since we would like to 
analyse naturally limited time series connected with two 
different phases of the solar cycle over the largest possible 
range of periods, we took the interval 43 to 1447 nHz (8- 
270 days) as the spectral window for all the investigated 
time series. The shortest period, 8 days, is connected with 
a possible range of correlations of the analysed data (see 
Table 2 and 3 where the autocorrelation coeficients with 
the lag = 1 day and lag = 7 days are given). The longest 
periods, 270 days, results from the actual length of the 
minimum time series equal to 540 days. The number of 
totally independent frequencies inside the chosen window 
is given by the value of the independent Fourier spacing, 
^/ifs = 1/?', where T is the time span of the data (Scargle 
1982). In the case of our time series, Tm — 540 days and 



ISN time series in the minimum phase. First, we calculate 
the Scargle power normalized by the variance of the data, 



A/ifsM = 21.433 nHz for the minimum, while for the rising 
phase Tr = 708 days and A/ifgR = 16.348 nHz. However, 
de Jager (1987) has shown by Monte-Oar lo simulations 
that the Fourier powers taken at intervals of one-third of 
the independent Fourier spacing are still statistically in- 
dependent. Thus, we accepted the numbers 198 and 259 
as the numbers of independent frequencies in the chosen 
window for the minimum and the rising phase time series, 
respectively. To illustrate the method used for determina- 
tion of the normalization factor k, we will process one of 
our time series step by step; as an example we take the 



(Horne & Baliunas 1986) for all 198 independent frequen- 
cies which allows us to construct the graph presented in 
Figure 2a. This shows the cumulative number of frequen- 
cies for which the Scargle power exceeds a certain value 
z. Then we fit all values of power z < 5 to the equation 
y = 198 exp(— z/fc), which gives the value of the normal- 
ization factor k equal 2.59, since when exp(— z/fc) << 1 
the formula for FAP reduces to FAP = A^exp(— z/fc). 
Finally, we normalize the power spectrum once more by 
dividing the Scargle power by 2.59 to obtain the normal- 
ized periodogram for which FAP values are easily calcu- 
lated. Therefore, if we substitute z = 30.54, fc ~ 2.59, 
and N = 198 in the FAP formula we get the normal- 
ized power NP = z^ = 11.79 and FAP = 0.0015 for 
the highest peak in the ISN minimum time series peri- 
odogram. Figure 2b presents this normalized periodogram 
of the minmum ISN series together with FAP signifi- 
cance levels obtained from FAP formula with = 198 
and fc 2.59. 

Since peaks in a periodogram may arise from aliasing 
or other phenomena not present in Gaussian noise (e.g., 
spectral leakage arising from the spacing of the data and 
from the finite length of the time series), the FAP val- 
ues alone, are insufficient for establishing whether or not 
strong peaks in a periodogram are indeed real periodici- 
ties in the time series. Also, some small peaks present in 
the original periodogram can be real in the case when the 
normalized factor fc would be too large in consequence of 
treating real periods as noise. We test for genuine peaks by 
recomputing the periodogram after randomising the data 
on the time g rid. T his procedure (Delache et al 



1985 



Ozgiig & Atag 1994 ) maintains the noise characteristic of 
the time series but destroys all coherent signals, especially 
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Fig. 3. a). The normalised periodogram of the minimum, original SVC 810/0 time series, b), c), ... same as a) but 
recalculated after successively removing from the original data one, two, and more sine curves having periods with 
peaks whose FAP values are smaller than 0.5%. In each graph the removed periods are indicated at the top. The 
dashed lines show FAP significance levels. 



those with periods longer than the chosen cut interval of 
the data. In our time series, we cut the data with a seven 
day interval. It preservs to some extent the correlation 
characteristics of the data, so if a period results from a 
strong correlation inside the data, the number of cases in 
which it is observed should be rather large. We repeated 
this simulation 10000 times, every time computing the 
number of cases in which the recalculated power values 
for the periods having peaks in the original spectrum are 



equal to or larger than the peaks power of the real data. 
The results of these calculations are presented in Table 4 
and 5 as small letters situated after the frequency of the 
peaks found in the original periodograms. The letter 'a' 
indicates that, for 10000 simulations, in less then 10 cases 
the peak value at the given frequency exceeded the corre- 
sponding peak power of the real data. Successive letters 
mark the intervals for which the probability (calculated 
from 10000 simulations) of obtaining as high a peak as in 
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Table 2. Charactaristic parameters for the minimum original time series(O) and all the next obtained after removing 
from the original data successively one, two and more sine curves ('-1', '-2', ...) having periods with peaks whose FAP 
values are smaller than 0.5 % These frequencies are printed in bold letters 
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Freq. 


Per. 
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var. 




[nHz] 


[day] 


1 


7 








cxp. 




[nHz] 


[day] 


1 


7 








cxp. 


1 


2 




3 




4 


5 


6 


7 


1 


2 




3 




4 


5 


6 


7 


ISN 


















MMF 





















84a 


138 


.839 


.126 


2.59 


11.8 


0.15 


17.4 





858a 


13.5 


.523 


187 


1.56 


19.0 


0.00 


12.6 




















-1 


419a 


27.6 


.476 


079 


1.54 


19.0 


0.00 


23.8 




















-2 


400a 


28.9 


.390 


095 


1.65 


13.7 


0.02 


30.0 


SVC 


405 
















RRE 


405 



















418a 


27.7 


.801 


.007 


2.50 


19.8 


0.00 


21.0 





416a 


27.8 


.714 


014 


2.25 


23.4 


0.00 


22.6 


SVC 


810 
















RR.E 


810 



















423a 


27.4 


.905 


.177 


1.56 


49.9 


0.00 


35.9 





418a 


27.7 


.801 


218 


1.46 


42.9 


0.00 


26.4 


-1 


396a 


29.2 


.895 


.285 


1.71 


18.7 


0.00 


45.2 


-1 


152a 


76.2 


.765 


346 


1.70 


27.7 


0.00 


40.6 


-2 


150b 


77.2 


.883 


.321 


1.89 


15.6 


0.00 


52.1 


-2 


381a 


30.4 


.707 


216 


2.06 


15.3 


0.00 


48.8 


-3 


76b 


152 


.874 


.228 


2.01 


14.1 


0.01 


57.7 


-3 


268b 


43.2 


.667 


242 


2.07 


11.8 


0.14 


54.2 


-4 


529b 


21.9 


.852 


.105 


2.10 


11.2 


0.28 


61.9 




















— 5 


238b 


48.6 


.836 


.164 


2.11 


11.5 


0.20 


66.0 




















SVC 


1215 
















RRE 


1215 



















424a 


27.3 


.912 


.219 


1.60 


34.6 


0.00 


26.7 





157a 


73.7 


.757 


271 


1.72 


18.8 


0.00 


13.7 


-1 


80b 


145 


.905 


.297 


1.72 


14.8 


0.01 


34.5 


-1 


420a 


27.6 


.718 


160 


1.91 


18.7 


0.00 


26.7 


-2 


151b 


76.6 


.893 


.207 


1.79 


15.0 


0.01 


41.8 


-2 


383b 


30.2 


.683 


221 


1.93 


12.2 


0.10 


34.3 


-3 


528b 


21.9 


.886 


.112 


1.90 


12.1 


0.11 


47.5 


-3 


267b 


43.4 


.649 


233 


1.89 


11.8 


0.14 


40.6 


-4 


113b 


102 


.873 


.176 


1.94 


11.5 


0.20 


52.9 




















-5 


498b 


23.2 


.855 


.084 


1.95 


11.8 


0.15 


57.8 




















-6 


242b 


47.8 


.839 


.129 


2.18 


11.8 


0.15 


62.5 




















SVC 


1620 
















RRE 


1620 



















425a 


27.2 


.898 


.111 


2.15 


23.1 


0.00 


25.0 





155b 


74.7 


.687 


197 


2.00 


11.7 


0.16 


10.6 




















^1 


420b 


27.6 


.653 


123 


2.07 


11.2 


0.28 


19.3 




















-2 


382b 


30.3 


.627 


159 


2.12 


10.8 


0.40 


27.2 


SVC 


2800 
















RRE 


2800 



















85b 


136 


.931 


.134 


2.38 


13.0 


0.05 


17.4 





383b 


30.2 


.647 


044 


2.32 


7.8 


7.78 


7.6 


-1 


424a 


27.3 


.921 


.020 


2.57 


13.5 


0.03 


28.0 




















SVC 


4995 
















RRE 


4995 



















90a 


129 


.862 


.286 


2.00 


18.3 


0.00 


14.3 





53a 


218 


.710 


347 


1.64 


23.5 


0.00 


19.4 


-1 


124a 


93.3 


.832 


.166 


2.20 


13.9 


0.02 


24.6 


-1 


158a 


73.2 


.625 


224 


1.83 


19.5 


0.00 


30.3 


SVC 


8800 
















RRE 


8800 



















122b 


94.9 


.635 


.355 


1.63 


15.8 


0.00 


10.3 





156b 


74.2 


.638 


396 


1.52 


16.2 


0.00 


10.0 


-1 


154b 


75.2 


.593 


.289 


1.71 


14.9 


0.01 


19.9 


-1 


121b 


95.6 


.598 


339 


1.54 


14.3 


0.01 


18.1 




















-2 


83b 


139 


.558 


285 


1.55 


11.8 


0.15 


24.3 



The columns of the table show: 1: time series('-l' indicates that the one sine curve with period given a row above was removed 
from the original data, '-2' indicates that the two sine curves with periods given in two rows above were removed from 
the original data, and so on), 2: frequency and period of the highest peak in a given time series, a small letter following a 
frequency value indicates to which interval of probability ('a': < 0.1%, 'b': 0.1 - 1%, 'c': 1 -5% ) this period belongs to after 
the randomising procedure, 3: the autocorrelation coefficients calculated with two lags equal 1 and 7 days respectively, 4: the 
normalization factor k computed according to the procedure described in Sect. 3, 5: the normalized power equal to the Scargle 
power divided by k, 6: the FAP value resulting from the normalized power, 7: in successive rows a percentage of the original 
data variance explained by the prominent sinusoidal signals found in a given type of data. 



the original periodogram by chance are as follows: 'b': 0.1 
- 1%, 'c': 1-5%, 'd': 5 - 15%. 

When more than one periodic signal is present in 
the data, multiple significant peaks appear in the peri- 
odogram. Alternatively, a true signal at frequency ujo can 
cause peaks in the periodogram at frequencies other than 
Wo because of the finite length of the data and irregu- 
larities in the data spacing. A useful procedure for de- 
termining whether any additional peaks with significant 
false alarm probability are physically real is an iterative 



peak removal technique (Delache & Scherrer 1983; Horne 
& Baliunas 1986| ). The highest peak in the original data 
periodogram provides the frequency ujo corresponding to 
the strongest sinusoidal signal present in the data. Using 
the method of least squares the phase and amplitude of 



this sinusoidal signal are fitted from the original data. This 
allows subtraction of the best fitted sine curve from the 
time series and then recalculation of a new periodogram. 
This procedure is repeated as long as the FAP of the peaks 
connected with subtracted sine curves are smaller than 
0.5 %, producing the main periods in each of the analysed 
time series. To see the effect of removing these peaks from 
the original data, we present in Fig. 3 the six consecutive 
periodograms, calculated according to this procedure for 
the case of the minimum SVC 810 time series. We see 
that the original periodogram is dominated by a strong 
peak at frequency 423 nHz (27'?4 ), so the structure of the 
periodogram near this frequency is difficult to recognise. 
When the sine curve with this period is removed from the 
data, the new period at 29'?2 begins to be visible (Fig. 3b). 
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Table 3. Same as Table 2 but for the rising phase data. 
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Subtracting the next sine curves from the data we come 
to Fig. 3f which presents the periodogram of the last de- 
scendant time series having period whose peak value gives 
FAP smaller than 0.5%. Comparing all the periodograms 
seen in this figure, we can conclude that the period values 
indicated by peak positions do not change much during the 
removed procedure, which confirms the significance of the 
periods found. The results of the removal procedure are 
summarised in Table 2 for the minimum and in Table 3 for 



the rising phase. Both tables show for the original time se- 
ries ( marked by 0), and those with successively removed 
sine curves (denoted by -1, -2, -3....) such characteristic 
parameters as: frequency and period of the highest peak, 
the autocorrelation coefficients, the normalization factor 
k, the normalized power equal to the Scargle power di- 
vided by k, the FAP value resulting from the normalized 
power. 
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Fig. 4. The normalised periodograms of the minimum ISN and SVC time series after removing from the original 
data all the sine curves having periods with FAP values smaller than 0.5%. The removed periods are shown at the 
top of each graph. The dashed lines indicate FAP significance levels. 
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Fig. 5. Same as Fig. 4. but for the rising phase ISN and SVC time series after removing from the original data aU 
the sine curves having periods with FAP values smaller than 0.5 %. 
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Fig. 6. Same as Fig 4. but for the minimum MMF and RRE time series after removing from the original data all 
the sine curves having periods with FAP values smaller than 0.5%. 
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Fig. 7. Same as Fig. 4. but for the rising phase MMF and RRE time series after removing from the original data 
all the sine curves having periods with FAP values smaller than 0.5 %. 
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Table 4. The characteristic periodicities (Unes) observed in the minimum. 
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The first nine columns give frequencies of tlic lines found in the indicated type of data. The small letter following the frequency 
denotes the interval to which belongs the probability to obtain by chance as high a peak at this frequency as in the original 
periodogram. The intervals are as folows: 'a': < 0.1 % 'b': 0.1 - 1 %, 'c': 1 - 5 %, 'd': 5 - 15 %. The bold letters indicate frequencies 
of the large peaks having FAP values smaller than 0.5 %. 

Columns 10 and 11 give the mean frequency and mean period of the line. They are calculated only from the periods marked by 
letters 'a', 'b', 'c'. The number under the mean frequency is the difference between the highest and the lowest frequency from 
all the periods found. 

The next three columns (12, 13, 14) describe in a shorter way in what type of data (magnetic, radio, spot numbers) the given 
line is visible. We mark this using the letters 'a', 'b', 'c', and two symbols '+' and The symbol '-' means that the given line 

is not noticed in the respective typo of data. The letters 'a', 'b', and 'c' indicate that the line is recognised and the periods 
found have the probability level described by these letters, while the symbol '+' indicates that in the given type of data are 
periods having the probability level described by the letter 'd'. The second row in column 14 gives the importance number 
which indirectly measure how large is support for this line in all the analysed time series (see text). 

The last two columns in the first row, determines the time series in which the highest peak for the given line have been observed. 
The second row gives the normalized power of this peak as well as the resulting FAP value. 



4. Results and Discussion 

Most of the calculated periodograms have been obtained 

from radio observations {SVC and RRE time series cal- 
culated for seven different frequencies), but the ISN and 
MMF data were also analysed. The final periodograms 
obtained when all the main sinusoidal signals given in Ta- 
bles 2 and 3 are subtracted from the original data are pre- 
sented in Fig. 4, 5, 6 and 7. Careful examination of these 
periodograms shows that many periods having a formal 
significance level near 30 % are seen in time series com- 



ing from different observational data. This fact allow us 
to suggest that even those periods could be real. 

Taking into account all the information about the pe- 
riods found obtained with the three different approaches 
to the data described in Sect. 3, we prepare Table 4 for 
the minimum, and Table 5 for the rising phase, which 
bring together information about characteristic periodic- 
ities (lines) observed in these phases of the solar cycle 
23 (we use the spectroscopic term 'line' for the mean pe- 
riod calculated from the periods observed in various time 
series, which, as we suppose, represent the same charac- 



S. Zi§ba et al.: Periodicities in solar data 
Table 5. Same as Table 4 but for the rising phase data. 
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teristic periodicity). As two time series were created from 
the radio data at each frequency, we use two successive 
rows in the radio data columns to separate periods ob- 
served in SVC and RRE radio time series. The exact fre- 
quencies of the detected periods were found with a 1 nHz 
resolution in the vicinity of peaks seen in periodograms 
constructed from normalized power values on the grid of 
the independent frequencies appropriate to the given win- 
dow. In all the figures, the points marked in the presented 
periodograms are situated in the places resulting from the 
grid of the independent frequencies, while the solid lines 
are drawn from the power values computed every 1 nHz. 
To discribe in some measure a 'strength' of the line we 
introduced so called the importance number of the line. 
This number shows how large is support for this line from 
all the analysed time series. It is calculated according to 
the following rule. The each letter 'a' in the line descrip- 
tion columns (columns 12, 13, 14) gives for the importance 



number '5', the letter 'b' gives '3' , and the letter 'c' only 
'1' . It is only one line with the importance number '5' in 
Tables 4 and 5. There are a few lines with a smaller im- 
portance number than '5' recognised in the analysed data, 
but we do not include them in Table 4 and 5. However, 
it is important to notice that some of the lines included 
in Tables 4 and 5 can be unreal in a sense that they are 
created from periods which in fact belong to two differ- 
ent but neighboring lines. For such a case the calculated 
mean period is somewhere between the periods of these 
two neighboring, unknown lines. The probability of such 
a situation increses for lines having long periods and large 
values of A/ (Column 10, second row). 

To aid in further discussion of the lines we have pre- 
pared Fig. 8, which shows all the lines present in Tables 
4 and 5. The level of darkness in this 3D graph illustrate 
the importance number of the lines. A close look at Fig. 
8 shows a clear difference among the 18 lines observed in 
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the minimum and the 22 hnes found in the rising phase 
of solar sunspot activity cycle. Although the 9 lines have 
almost the same periods in both the phases, the strength 
(measured indirectly by the importance number) of nearly 
all the 18 lines change, indicating that perhaps the physi- 
cal mechanisms responsible for them also change with the 
solar phase. From Fig. 8 it is evident that the lines in the 
rising phase gather into three groups: 



1. 



Magnetic lines: all lines with periods shorter than 16 
days. They have small importance numbers and were 
included on our list only because of very high peaks ob- 
served at these short periods in the MMF (mean mag- 
netic field) time series. In the minimum, only one line 
(13^5) from this group is enough strong to be present 
in Fig. 8, but in the minimum MMF periodogram (see 
Fig. 6) some of them are easily seen. 
Rotational lines: all lines with periods inside the 25 - 34 
days time interval of the Sun's rotation. The three lines 
from this group are visible in the minimum. The line 
27^3 has the largest importance number and results 
probably from the rotation of such solar phenomena 
as new solar cycle sunspots as well as the long-lived 
coronal streamer structure observed during the 1996 
minimum (Lewis et al. 1999). In the rising phase the 
period of this line shifts to 27'?7. The next two lines 
29'? 1 and 30?3 are stronger in the minimum than in 
the rising phase. This additionally supports the suppo- 
sition that they are associated with some medium, and 
large scale magnetic structures (coronal neutral sheet, 
global neutral lines) which dominate in the minimum 
and then decline (Lantos et al. 1992; Lantos 19991 ). 



In the rising phase four new lines appear (25.2, 26r0, 
26?8, 31?8). Three have periods shorter than those dis- 
cussed above. One of them, having a period of 26^0, is 
the third dominant line of this phase. We suppose that 
this line is strongly connected with eruption of new 
active-region magnetic fields emerging within the com- 
plexes of activity ' activ e longitudes' (Bumba & Howard 
|l969 ; Ruzmaikin 1998| ) and its period may be identi- 
fied with a fundamental period of unknow Sun's clock 
which value is freqcuntly taken as equal to 25'?5 or 25'?8 
(Bai & Sturrock |1993| ; Bai |1994D . 
3. Activity lines: all lines with periods longer than 75 
days. We propose this name as the strongest 151 days 
line in this group is thought to be related to the well 
known periodicity near 154 days seen mainly in the 
occurrence of high energy solar flares at the maximum 
of the solar sunspot activity. It is interesting to notice 
that this line is also observed in the minimum, but 
with a small importance number and only in the radio 
time series. 

In the minimum, the distribution of the line periods in the 
investigated time window is much more uniform. There are 
5 lines with periods shorter than 25 days but only one be- 
longs to magnetic lines. Within an accuracy of ±0?3 all 



of them were noted before in various time series (Hughes 
&:Kesteven [l98lt Pap et al. |1990D . The lines with mean pe- 
riods 13'?6 and 23^2 are the most interesting. They belong 
to the group of 9 lines which are observed in both phases 
of the solar cycle investigated. Although in our data their 
importance numbers are not large, they were recognised 
as the dominant periods in the SMM/ACRIM total irra- 
diance data for the years 1980 - 1988 (see Papet al. p^90| ). 
A very intriguing group of lines is formed by 4 lines (39, 
43, 48, 58 days) also observed in the minimum. All these 
lines are seen only in the time series which were obtained 
from the radio data, except for one with period 58 days 
also observed in the ISN data. Among them only the line 
at 39 days is also present in the rising phase, although a 
small shift of its period is visible. Coming to longer peri- 
ods we see that most of them are present in both investi- 
gated phase of the solar cycle. However, in the rising phase 
we observed more lines and all of them are stronger than 
those in the minimum. In the minimum, the strongest line 
in this group has period 75 days. However, in the rising 
phase this line is more weaker. It is not clear if it should 
be correlated with periods 73-78 days observed by some 



autho rs (Bai & Sturrock |l99l| ; Bai |19924 Ozgiig & Atag 
1994) mainly in flare activity during the maximum of a 
solar cycle. We suppose that this line could be important 
if it is a real harmonic of the best known line with pe- 
riod near 154 days. It is important to notice that all of 
them were mentioned before in various analyses (Lean & 
Brueckner Il989| Pap et al. |1990|; Bai & Sturrock |1991 



Kile & Oliver |1991| ; Antalova ^l999| ) 



Many previous studies by a number of authors have 
resulted in a wide range of solar periodicities, which are 
not easy to explain. This indicates that the problem of 
solar periodicities is still open and more systematic efforts 
should be undertaken. Here, we do not want to discuss 
all possible causes of the observed periods, but we want 
to present a suggestion which may be of help in further 
investigations. 

Recently Oliver et al. ( 1998D proposed that the peri- 
odic emergence of magnetic flux, manifested as sunspots, 
triggers the near 158 day periodicity in high-energy solar 
flares. As different magn etic features have di fferen t rates of 
rotation (Oilman |l974 ; van Tend & Zwaan 1976 Erofeev 
1999) we think that a periodic emergence and a constant 
conversion of various magnetic structures explain the ori- 
gin of the observed lines and their transformation with the 
phase of solar cycle 23. The main arguments supporting 
this idea are as follows: 

1 . There is a clear difference between the lines obtained in 
the minimum and in the rising phase. The lines char- 
acterictic for the minimum are probaly connected with 
both the short-lived small-scale magnetic fields which 
ori ginate fairly high in the convective zone (Oolub et 



al. |1981| ; Rivin |1999[ ) and large-scale fields structures 
like coronal holes, coronal neutral sheet, global neutral 
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Fig. 8. Characteristic lines generated by solar phenomena 
in a) the minimum and b) the rising phase of solar cycle 23. 
The three group of lines are inicated: The level of darkness 
qualitatively illustrate the importance of lines. 

lines. All these small and large scale magnetic fields 
dominate in the minimum phase of solar cycle (Golub 
et al. |1981| ; Wang et al. |l99^ ; Lantos [19991 ). 

2. In the rising phase a new distinct rotational line at 
26.0 days is seen. This period is equal to the rota- 
tion rate of the active longitudes (zones) at 30° (where 
the new flux appeared) determined from SOHO/MDI 
magnetograms by Benevolenskaya et al. ( |1999| ). They 
found in the northeren hemisphere the rotation rate 
446.6 ± 1.7 nHz and 444.8 ± 1.6 nHz in the southern 
one, while our value for the whole Sun data is 444.4 
± 4 nHz. Also the rotation rate of the equatorial zons 
equal 461.8 nHz from SOHO observations is almost the 
same as the period 25'?2 (459.5 nHz) seen in our radio 
data (see Table 5). 

3. The line at 151 days is in the rising phase very p romi- 
nent. This period is connected by Bogart ( 1982 ), Bai 
(|1987|) and Oliver et al. (|l998[ ) with strong, long-lived 
active regions giving most of the energetic solar events. 
This regions are generated by the emergence within the 
'active longitudes' magnetic tubes formed in the dy- 
namo located at the base of the c onvec tion z one (G olub 
et al. |l98l| ; Benevolenskaya et al. |l999| ; Rivin |l999D . We 
also observe in the rising phase ISN periodogram pe- 
riods such as 13'?lor 23^8, which were linked with the 
'active longitudes' and active centers in deep layers by 
Bai (119871). 



5. Conclusion 

To summarise, our analysis of radio observations gives the 
majority of the known periods, reveals a clear difference 
among periodicities observed in the two phase of the so- 
lar cycle examined and shows that daily measured radio 
fluxes at various frequencies are very useful for the system- 
atic study of solar periodicities observed in the different 



layers of the solar atmosphere. We are preparing a similar 
analysis for the next phases of solar cycle 23: the maxi- 
mum and the declining phase. 

We found the rotation rate of the 'active longitudes' 
in the rising phase as equal to 444.4 ± 4 nHz (26^0 ± 
0'?3). We suppose that this period can be identified with a 
fundamental period of unknow Sun's clock as a lot of the 
known periodicities are subharmonics of it. We think that 
in the minimum the lines observed are conected with the 
small and large-scale magnetic fields which then dominate 
, while in the rising phase most lines are generated by 
new magnetic structures connected with long lived active 
regions formed within 'active longitudes'. 

To understand the cause of all the observed periodici- 
ties, we still need new observational data, but more careful 
analyse of old radio data can also be useful (in prepara- 
tion). However, our investigation, idicates that the solu- 
tion of the observed solar periodicities should be sought 
in a complicated Sun's magnetic system which generate in 
the different solar data the compound set of solar period- 
icities. 
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